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Abstract 

We establish rigorous lower bounds on the speed of traveling fronts and on the bulk burning 
rate in reaction-diffusion equation with passive advection. The non-linearity is assumed to be of 
either KPP or ignition type. We consider two main classes of flows. Percolating flows, which are 
characterized by the presence of long tubes of streamlines mixing hot and cold material, lead to 
strong speed-up of burning which is linear in the amplitude of the flow, U. On the other hand the 
cellular flows, which have closed streamlines, are shown to produce weaker increase in reaction. 
For such flows we get a lower bound which grows as U 1 ^ for a large amplitude of the flow. 

1 Introduction 

Propagation of thin fronts in moving fluids arises in many situations in physics and engineering. 
Consider a mixture of reactants interacting in a region that may have a rather complicated spatial 
structure but is thin across. The reaction front moves towards the unburned reactants leaving behind 
the burned ones. When the reactants are mixed by an ambient fluid then the burning rate may be 
enhanced. The physical reason for this observed speed-up is believed to be that fluid advection tends 
to increase the area available for reaction. Many important engineering applications of combustion 
operate in the presence of turbulent advection, and therefore the influence of advection on burning 
has been studied extensively by physicists, engineers and mathematicians. In the physical literature 
one can find a number of models and approaches that yield different predictions - relations between 



the turbulent intensity and the burning rate pL 19, 2C, pj]. These results are usually obtained using 



heuristic models and physical reasoning. For a recent review of some of the physics literature we 



refer to |27|, §9|. 



The lack of agreement between different physical models makes rigorous results, even for simpli- 
fied mathematical models, particularly valuable and useful. A well-established mathematical model 
that describes a chemical reaction in a fluid is a system of two equations for concentration C and 
temperature T of the form 

v 2 

T t + u-VT = KAT + ^g(T)C (1) 

K 

2 

C t + u-VC = ^AC - %(T)C. 
Le k 

For exposition purposes, all consideration in this paper will be carried out in two spacial dimensions, 
but our methods extend to an arbitrary dimension in a straightforward way. Equations (||) are 
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coupled to the reactive Euler equations for the advection velocity u(x,y,t). Two assumptions are 
usually made to simplify the problem: first, constant density approximation || that allows to 
decouple the Euler equations from the system (|l]). Then one may consider u(x,y,t) as a prescribed 
quantity that does not depend on T and C . Furthermore, it is often assumed that Le = 1, or, 
equivalently, thermal and material diffusivities are equal. These two assumptions allow to reduce 
the above system to a single scalar equation for the temperature T: 

r)T v 2 

— + u(x, y, t)-VT = kAT + ^/(T) (2) 

with f(T) = g(T)(l — T), provided that C(x,y,0) = 1 — T(x,y,0). We will consider the problem 
(H) in a strip Q = M x x [0, H] y with boundary conditions in x: 

T(x,y,t) — > 1 as x — > — oo, T(x,y,t) — > as x — > +oo (3) 

and either Neumann 

dT dT 

%-(x,0,t) = %-(x,H,t) = (4) 
dy dy 



or periodic 



T(x,y,t)=T(x,y + H,t) (5) 



boundary conditions in y. Furthermore, we assume that the initial data To(x,y) for (0) satisfies the 
bounds 

T (x,y) = 1 - 0(e Xx ) for x < 0, T (x, y) = 0(e" Az ) for x > 0, (6) 
jVT | = 0{e- x \' x \) for some A > 0. (7) 

We adapt fairly general assumptions on /, requiring only that f(T) is not equal identically to 
zero, and 

/(0) = /(l)=0, f(T)>0 for TG (0,1), f € C^l]. (8) 

Two types of reaction rates f(T) are distinguished in this class. The KPP-type reactions satisfy 

/(0) = /(l) = 0, /(T) > for T £ (0, 1). (9) 

An additional requirement f'(0) = maxygfo i]/(T)/T is often made. We do not make such require- 
ment in this paper, and call the class described by (||) general KPP. Our interpretation of KPP 
includes an important Arrhenius-type non-linearity, 

f{T) = C{l-T)e- A ' T , 

that is believed to be an appropriate model for many chemical reactions in the context of reaction- 
diffusion models. We also consider the ignition non-linearities with 

f(T) = for T € [0, 9 ] and T = 1, f(T) > for T £ (0 , 1). (10) 

By our assumptions on the nonlinearity (|8|), we can find 9 4 > 9\, and /o, C > such that 

f(6)>f far 6 (0i -C, 04 + 0- (11) 

The values of the constants /o, C an d $1,4 are the only information on the nonlinearity f{T) that 
shows up in our bounds on the burning rate. 
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We assume that advection u(x, y) G C 1 (f2) is time independent, has mean zero in the x-direction: 

H 

Ul {x,y)dy = (12) 



o 

and is incompressible: 

V-u = 0. (13) 

The mathematical literature on the scalar reaction-diffusion equation @ is enormous; far from 
giving an exhaustive overview, we mention several papers directly related to our work. First rigorous 
results about traveling waves for equation @ go back to classical works of Kolmogorov, Petrovskii 
and Piskunov |^3| and Fisher |13||, which considered the case u = in one dimension for the KPP 
nonlinearity. Recently equation (|2|) with and in particular the effect of advection, became a 

subject of intense research. Berestycki and Nirenberg [||, (7|, and Berestycki, Larrouturou and Lions 
0] initiated the studies of the existence of traveling waves for equation (||) of the form 

T(x,y,t)=T(x-ct,y), (14) 

for shear flows of the form u = (u(y),0). Their stability was studied in ||, 25, pCf| , while in further 
works [||, [}l], 32] stability and existence of traveling waves were established for the wider class of 
periodic flows. In this case, the traveling fronts have the form 

T{x,y,t) = U(x - ct,x,y) 

and are periodic in the last two variables. These and other results were recently reviewed in f33[1 , 
and we refer the reader to this paper for a detailed exposition of the subject. Until very recently, 
there were no rigorous results on the physically interesting question of the speed of traveling waves. 
First such results have been established in || for percolating flows, and in (l| and [^] for the shear 
flows. Numerical studies of the propagation of fronts were performed for a shear flow in [21| with 



Le 7^ 1, and for cellular flows in [22]. 



Another major direction of research has been homogenization approach. The homogenization 
regime k — > 0, when the front width goes to zero, was extensively studied for KPP-type nonlinearity 
and for advection velocity that is periodic and varies either on the integral or diffusive scale by 



Freidlin []14j, |l6fl . Recently Majda and Souganidis derived an effective Hamilton- Jacobi equation 



in the limit k — > for the case of advection velocity varying on a small k— dependent scale that 



is larger or comparable to that of the front width 24]. This effective equation is still difficult to 



analyze, and analytical predictions have been derived only for the shear flows. Numeric experiments 



exploring the results of |24] h ave been carried out in yl], 11, pq ]. 

Very recently, Hamel |17| and Heinze, Stevens and Papanicolaou |l8| proposed an elegant varia- 
tional approach to the estimates of the speed of traveling waves in the presence of periodic advection. 
However, to the best of our knowledge nontrivial lower bounds using this method were obtained so 
far only for shear flows in the homogenization regime or for small advection, where they provide 
precise bounds for the small speed-up of the front 

The key question we wish to address in this paper is: what characteristics of the ambient fluid flow 
are responsible for burning rate enhancement? The question needs first to be made precise, because 
the reaction region may be complicated and, in general, may move with an ill-defined velocity, when 
traveling fronts do not exist. To measure the speed of burning in such situations, the bulk burning 
rate 

V(t) = J T t {x,y,t) d -^- (15) 
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Figure 1: The structure of the shear flow 



and its time average 

t 

(V)t = \J V(s)ds 
o 

have been recently introduced in Q. Note that for traveling fronts of the form ( |i~4| ) we have V(t) = c, 
but the notion of bulk burning rate makes sense in much more general situations when traveling 
fronts of the above form may not exist, and bulk burning rate serves as a natural generalization of 
the front speed. We have obtained in || lower bounds for (V)t when f(T) is a concave function of the 
KPP type. The bounds are linear in the magnitude of the advecting velocity u(x,y) provided that 
there exist tubes of streamlines that connect x = — oo and x = +00, satisfying some mild additional 
technical assumptions. We say then that the flow is percolating. In particular these bounds hold for 
shear flows of the form (u(y),0). 

In this paper we consider much more general reaction rates f(T) that are either of the ignition or 
general KPP type, and establish similar lower bounds for V(t) for percolating flows that are periodic 
in space. The bound is linear in the magnitude of u and deteriorates as the scale of oscillations of u 
becomes comparable to the laminar front width I = k/vq. It is easy to show ||] that for any it G C , 
the burning rate (V)t satisfies linear in ||u||oo upper bound (for initial data as in @ and (0)). 
Therefore, shear (and, more generally, percolating) flows are as effective as possible in speeding up 
combustion in terms of the power of || lijloQ in the large intensity regime. In particular, we establish 
the following lower bound for the bulk burning rate in a shear flow. 

Theorem 1 Let To(x,y) be an arbitrary initial data satisfying l^) and ffi), and let T(x,y,t) satisfy 
^) with the either the Neumann fy) or periodic boundary conditions (^). Let also u(x,y) = (u(y),0) 
in Then both for KPP and ignition non-linearities we have 

lim<n>c(,„ + £(l + i-) 1 / (16) 
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Figure 2: Streamlines of a cellular flow {H = ir). 



where the constant C depends only on the reaction function f . Here the intervals Ij = [cj — hj , Cj + 
hj] € [0, if] are any intervals such that 

\\ u \\oo,J < \u(y)\ < \\u\\oa,j, \\u\\aa,3 = SUp|u(y)|. (17) 

2 h 
We do not require Uff,- = [0, If]. 

The choice of intervals Ij is up to us, and should be made to maximize the lower bound. See Figure |] 
for an illustration. 

As a corollary, the bound (|l6| ) holds for the speed c of a traveling front of the from (|l4|). Our 
bound behaves correctly in the homogenization regime when u(y) has the form u(y) = j-v(yfe) and 
provides a bound that is linear in the magnitude A of advection, in agreement with || [18|], where 
homogenization limit was studied. We also prove the analog of Theorem [l] for general percolating 
flows (see Theorem || in Section ||) . 

Another main result of this work concerns cellular flows with closed streamlines. Roughly speak- 
ing, in terms of their burning enhancement properties, such flows can be thought of as "the worst" 
class of flows, opposing "the best" percolating flows. One can expect the burning enhancement to 
be significantly weaker for cellular flows because of the numerous diffusive interfaces which prevent 
hot and cold regions from mixing fast. Cellular flows pose mathematically more challenging problem 
because of these diffusive interfaces; we will see that the estimates for percolating flows will form 
only a fraction of the argument we will need in the cellular case. We consider a particular example 
of a cellular flow 

u( x , y ) = UH<^-,-^y ^ y ) =s i n (|) sin (|). (18) 



See Figure g for an illustration of streamlines in a single cell. Then under the assumptions of a large 
Peclet number and thin laminar front width: 

Pe = — >1, -<1, l = - (19) 

K H Vq 

we get the following result. 

Theorem 2 Let To(x,y) be an arbitrary initial data satisfying ^) and ffi), and let T(x,y,t) satisfy 
(fl) with the either the Neumann fy) or periodic boundary conditions (^). Let also u(x,y) in (11) be 
given by fiTl\), and assume is satisfied. Then both for KPP and ignition non-linearities we have 

\ (Cuf^ + C 2 )v , if T C <T U 
lim (V) t >{ / / N 1/5 \ (20) 

t— +00 II.., . 

Here r c = k/vq is the chemical reaction time and t u = H/U is the turnover time. The constants in 
the inequalities depend only on the reaction f , more particularly on the constants /o, C) an d $4 — $i 




that appear in 



To the best of our knowledge this is the first rigorous bound on the traveling front speed in a cellular 
flow. Note that the change of behavior of our bound depending on the ratio t c /t u is physically 
natural since for t u <C t c the front folds onto itself inside the period cell, which diminishes the affect 
of advection. The lower bound of Theorem ^ displays square root dependence on the flow intensity 
U until U reaches a critical value determined by a condition r c = r u . After that, the lower bound 
behaves like U 1 ^ . Recently, Audoly, Berestycki and Pomeau [Q gave an heuristic argument which 
proposes that the speed of the traveling front for cellular flows should scale as U l / A in the large U 
limit, which may indicate that our lower bound is not far off from the sharp bound. 

One of the fundamental mathematical difficulties we deal with in this paper may be roughly 
described as follows. We will be able to bound the burning rate from below by integrals over the 
domain of reaction term f(T) and of the square of the gradient |VT| 2 . It turns out that in order to 
obtain a lower bound on V in terms of u, it will be necessary to bound the integral of the higher 
derivative Laplacian term in terms of integrals of f(T) and |VT| 2 . One can expect to do this using 
parabolic regularity, but the constants in such a priori estimates typically depend on u, and this 
dependence turns out to be too crude to get interesting results. We tackle this difficulty by taking 
advantage of the fact that what we need to estimate is the integral of Laplacian, not of the absolute 
value of Laplacian, and employ an appropriate averaging procedure to reduce derivatives. We hope 
that this idea will be useful in other related contexts in PDE estimates. 

The paper is organized as follows. We prove Theorem |l] in Section ||, as well as some other results 
for shear flows. The analogous results for the percolating flows are proved in Section 0. We present 
our main results for the cellular flows, in particular implying Theorem ^ in Section [|. Sections ^, 
^ and H contain some intermediate estimates in the proof of Theorem |2| We put these estimates 
together to finish the proof in Section ^. 
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2 Shear flows and general nonlinear it ies 



We first consider (|2|) in a shear, or unidirectional, flow (u(y), 0), which is a particular example of a 
percolating flow. The proofs are somewhat less technical in this case and allow us to introduce some 
of the ideas used in the general case. Equation @ in a shear flow becomes 

§ + u W f = .AT ,2!) 
T{x,y,0) = T (x,y). 

The advection is assumed to be mean-zero: 

n 

u{y)dy = 0. (22) 



We impose an additional assumption 



dT 

m a a (23) 



This condition is satisfied for all times provided that initially we have 



K ATo-u(y)^ + -^f(T )>0, (24) 



as follows from the maximum principle (see, e.g. |28|]). Therefore fl23|) is not a constraint on the 
dynamics but rather on the initial data. We note that ( p3| ) is true for traveling fronts of the form 
T(x — ct,y) both for general KPP and ignition nonlinearity f?j. We assume the usual boundary 
conditions @ at the left and right ends of the strip f2 and either Neumann (|j) or periodic (p|) 
boundary conditions at y = 0, H. We also require that the initial data Tq(x, y) satisfies (y) and (0). 
These conditions are preserved by evolution (see e.g. ||) if the advecting velocity u(y) £ C [0, H], 
that is, we have for each t > 0: 

1 - T{x,y,t) < C(t)e Xx for cc < 0, T(x,y,t) < C{t)e~ Xx for x > 0, (25) 
\VT(x,y,t)\<C(t)e- x W 

provided that (||) and (|7|) hold initially. 

Let Ij = (cj — hj,Cj + hj) C [0, H] be a collection of intervals satisfying (0). In particular u(y) 
does not change sign on the intervals Ij. We do not require that Djlj = [0,H]. Then the bulk 
burning rate V{t) defined by (|l~5|) obeys a lower bound described by the following Theorem, which 
is the first main result of this section. 

Theorem 3 Let T(x, y, t) be a solution of with the boundary conditions ^) and either fy) or 
([|). Let the initial dataTo(x,y) satisfy p^), ^) and Q). Furthermore, assume that u(y) £ C [0, H] 
has mean zero ^2^) and the nonlinearity f(T) satisfies ^). Then there exists a constant C > that 
depends on f(T) but not on To(x,y) oru(y), such that for any collection of intervals Ij that satisfies 
^P\ ) we have 

V(t)>C wo + Wl + lj / |u(y)||M (26) 

■ ii K 
with I = — . 

VQ 
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Remark. 1. The lower bound ( p6|) does not deteriorate when oscillations in u(y) become faster in 
space as long as its amplitude grows according to 



v 




1 



This agrees well with the homogenization limit u e (y) = —u(y/e) considered in || and [18], that 

produces speed-up of the front of order 0(vq). This is also an improvement of the analogous lower 
bound for V(t) for the convex KPP case obtained in [0], where (l/hj) 2 appeared in the factor. 

2. The regularity assumption on u(y) is used only to guarantee preservation of the boundary 
conditions (|7|) that allows us to integrate by parts in the proof. None of our bounds depend on the 
size of derivatives of u(y). 

It has been shown in [Q, |7| both in the case of ignition non-linearity (|l0|), and for the general KPP 



nonlinearity (]9D that there exist traveling front solutions of (21) of the form T(x,y,t) = U(x — ct,y). 
The speed c = c* is uniquely determined by the nonlinearity f(T) and advection u(y) in the ignition 
case, while traveling front solutions exist for c > c u for some minimal speed c u in the KPP case. 
The function U(s,y) is monotonically decreasing in the variable s = x — ct in both cases, so that 



(23) holds. Theorem g implies the following estimate on the speeds of the traveling fronts. 



Corollary 1 Let T(x — ct,y) be a traveling front solution of (21) with f(T) being either of the 
ignition nonlinearity type JT^J, or of the KPP type Then there exists a constant C > that 
depends on the function f but not on u(y) such that 



c>C 



1 ' 1 \ f i / m d v 



V 



-o + E|l + X-J J l«G0lf| (27) 



Corollary |l] follows immediately from Theorem [3| since we have V(t) = c for T(x, y, t) = U(x — ct, y) 
due to the boundary conditions U(s, y) — > 1 as s — > — oo, U(s, y) — > as s — > +oo. Corollary 1 and 
the stability results for traveling fronts ^fj, ^2| imply Theorem [l] for general initial data. We prove 
now Theorem |j] assuming the result of Corollary |]. 

Proof of Theorem |l]. Consider first ignition non-linearity. In this case we will show that 

lim (V) t = c*, 

t— >oo 



where c* is the unique speed of the traveling front. Then (16) will follow from (|27|). It was shown in 



[32 1 that for the initial data satisfying (g) (actually just tending to 1 and at the two ends) there 
exist functions £1,2 such that 

\£i(t)\ = o(t), as t -> 00, 
and functions qi(t,x) that satisfy the linearized problem 

% , ,dqi 

— + u(y)- = A qi 

such that 

U(x - c*t + fi(£),2/) -qi(x,y,t) <T(t,x,y) < U(x - c*t - &(t),y) + q 2 {x,y,t). (28) 
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Here U(x — c*t,y) is the traveling wave solution of fl2l|). The initial data qi(x,y,0) may be chosen 
in L 1 n L°°(0) Then we have for any c > c*: 

(V) T = i / dt J ^^-T t dxdy = ^J ^^-[T(x, y, r) - T (x, y)] 
on n 

OH cr H 

= \ J dxj ^[(l-T )-(l-T)] + i JdxJ ^L[T(x,y,T)-T (x,y)} 

-oo 

1 °° H d 

+- I dx J ^[ T 0>y> r ) - T o( x >y)] 

CT 

oo H oo H 

C 1 f , /" r . _ . . , 1 f 1 f dy . C 

<— + c+- / dx / — U(x - ^2{t),y) + - / dx / — q 2 {x + c*t,y, t) < — + c 

(c-c)r cr 

and hence limsup(V) T < c*. Similarly one may show that liminf(V) T > c' for any d < c*, which 

T^OO ' ' T^OO 

shows that (16) holds. 

In the KPP case, the estimates of the sort (28) are not yet available. However the bound (16) 
can be shown by reduction to the ignition non-linearity case. Indeed, given KPP type reaction /, 
consider ignition type reaction /e = Xe/ < /, say by cutting / off in a small neighborhood near zero 
for T < 9. The constant C in (^) does not depend on for 6 small enough as will be seen from the 
proof of Theorem ||. Let T and Tg satisfy equations with reactions / and fg respectively with the 
same initial data Tq(x, y). Then Z = T — Tg satisfies 

Zt + u ■ VZ - kAZ = V \f(T) - fg(T e )) > V l{f{T) - f(Tg)). 

K K 

It follows from the maximum principle that if Z(x, 0) > then Z(x, t) > 0. Hence for the same initial 
data, the burning rate for the KPP reaction / is not smaller than for the ignition non-linearity fg: 

(V[T]) t = \J ^[T(x,y,t)-T (x,y)] > -J ^-[T e (x, y, t) - T (x, y)] = (V[Tg]) t . 



This implies the validity of the lower bound ( |16| ) in the KPP case. □ 

We now turn to the proof of Theorem ||. The proof follows the general ideas of || with significant 
modifications required since / is not of concave KPP class. Our starting point is the following 
observation. 

Lemma 1 Under assumptions of Theorem [| we have 

V(t) = |/ f(T(x,y,t))^ > K J \VT(x,y,t)\^. (29) 



Proof. The equality in (^) is obtained simply by integrating ( |2lD over Q using the boundary 
conditions d3|) and (Q) or (j^) , and mean-zero condition ( p2[ ) on advection. To get the inequality we 
multiply ([n]) by T and integrate over Q to get 

jTrfSL + K j ivrfx^.t)!^ = i / r /<r)** < v(t). 
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This implies (gg) since T t > and < T < 1. □ 

As a warm-up, we now prove a simple and general proposition, which already provides a glimpse 
of some of the ideas which we will use to obtain more sophisticated results. Namely, we show that 
for any divergence-free velocity u(x, y) satisfying mild regularity conditions (it does not have to be 
a shear flow), and solution T satisfying (p3|), the burning rate is bounded below by Cvq. 



Proposition 1 Let T(x,y,t) be a solution of (\^) with the boundary conditions ffi) and either ftfy 
or t^). Assume that u(x,y) G C 1 ([0,ff] x M) satisfies (fZl j and and that non-linearity f(T) 

satisfies ^). Let the initial data To(x,y) satisfy (2£), $Q)and (Qj. Then there exists a constant C, 
depending only on the parameters £ and /o, such that 

V(t)>C((,f )v 

with the constants ( and /o defined in (|7"/|j. 

Proof. The proof is similar to the proof of Lemma 2 in Q . We can find y such that 

H 

J \VT(x,y)\ 2 dx < ^ J J \VT(x,y')\ 2 dxdy', 



and 



M 



H 



J f{T[x,y))dx <^Jj f(T(x,y'))dxdy'. 
m o m 

Then we can find x\ < X2 such that T(x±,y) = 64 + £, T(x2, y) = 9% — T(x, y) G [^1 — C> #4 + C] if 
x\ < x < X2 (see (111]) for the definition of #1,4). Then 



|VT(x,y)| 2 dx > 



JR 



h~0i + 2C 

|X2 - Xl| 



and 



Therefore we have 



Hence we obtain 



f( T (x,y))dx > / |x 2 - xi| 



\ 



H 



H 



J J \VT\*dxdy J J f(T)dxdy>Cfl l2 QH. 



M 



H 



M 



M 



k|VT| 2 + —f(T) 

K 



—jj - > Cf ' (v . 



Then Proposition |^ follows from Lemma |l|. □ 

We now return to the shear flows. To obtain more precise bounds involving advection velocity 
u(y), we will bound from below in terms of u(y) either the integral of f(T) or the L 2 -norm of |VT|, 
and use Lemma |l|. The general plan in || was to integrate over all axis in x, obtaining an equation 
with an explicit term u(y) in it. We were able to bound the rest of the terms from above by a 
combination of / f(T) and J | VT| 2 after averaging in y and t to control Tt and AT. 

An additional twist we need here is to reduce our consideration to the region in space where the 
reaction actually takes place. In the case of ignition non-linearity, there is no reaction for sufficiently 
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low temperatures. Similarly, for the Arrhenius type non-linearity, reaction is very weak near T = 0. 
On the technical side, restriction of consideration to some finite time dependent domain D with 
T in appropriate range will mandate additional averaging in x to control all terms by J f(T) and 
/ |VT| 2 . We will identify a region D in x such that on one hand the temperature has a certain drop 
over this region and on the other for every x E D there is some y E Ij such that reaction is bounded 
away from zero at the point (x, y). This will provide us with two alternatives for each x E D: either 
reaction is uniformly bounded away from zero for that x or temperature drops by a certain amount 
on the interval x x Ij. In the first case J T f(T)dy will have to be large and in the second fj \T y \ 2 dy 
will be bounded from below. Then we will integrate (|2l]) over x E D at a fixed time t. That will 
relate u(y) to some terms involving V(t) and AT. We will additionally average both in x and y, 
which will bring AT into a form that can be bound by a combination of integrals of /(T) and | VT[ 2 . 
That will be possible since these have to be large on D as explained above. Finally we will use 
Lemma |l| to finish the proof. 

In order to define the region D where much of the reaction takes place let us fix #4 > #3 > #2 > Q\ , 
where #4, 9\ are as in (|TT|) . Let Ij be an interval on which ([l?]) holds with u(y) > (the case of Ij 
where u(y) < is similar). We fix time t > and choose two points xo an d x\\ 

xq = inf {x : for any x' > x there exists y E Ij s.t. T(x', y, t) < 64} 

x\ = sup {x : x > xo and for any x' E (xo, x) there exists y E Ij s.t. T(x', y, t) > 9±} . 

In other words, for any x E [xo,xi] = D there exists y E Ij such that T(x,y,t) E [^1,^4], and hence 
f(T(x,y)) > fo- Note that xo is well-defined and finite since T(x,y,t) — > 0, 1 as x — > ±00 uniformly 
in x because of (|25|). The definition of xq implies that 



T(x ,y)>8 4: for ally (30) 

and thus x\ is well-defined. Moreover, 

T(x u y) < 9i for ally. (31) 

In preparation for multiple averaging in y that will be performed to control AT let us introduce 
the function G(h,£): 

I 0, iei > h 

that corresponds to the following averaging in y: 



dS / p(y)dy= / G(hj ,y - Cj)p(y) dy (32) 



1 

hj 

for a test function p(y). The two integrations when applied to AT are required to get rid of derivatives 
of T. Observe that the function G(h, £) has the following properties 

0<G(fc,0<l, G(h,0>l for (33) 



We note that (pQ) and (pl|) imply that at the two ends of the interval [xo,xi] we have 

3 

dyG(hj,y - Cj)u(y)T(x ,y) > 6 A Fj 
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Cj+hj 



dyG(hj,y - c j )u(y)T(x 1 ,y) < 6 1 F j 



with 



Cj+hj 



dyG(hj,y- Cj)u(y). 



In preparation for averaging in x, we choose rjo and rji so that 

{Cj+hj 
£ > : J dyG(hj,y - Cj)u(y)T(x Q + £,y) = 8 3 Fj 
Cj-hj 

! Cj+hj 
£>0: J dyG(h j: y-c j )u(y)T(x 1 -Z,y) = 2 F j \ . 
c j- h j ) 
We remark that xq < xq + tjq < x\ — rji < x\. Now we are ready to average (El). Given a G (0, 770 



and j3 E (0, rji) we integrate (21) in x E (xq + a, x\ — (3) and in y according to (32]): 

XI— P c 3+ h j XI- P c j+hj 

J dx J dyG(hj,y - Cj)T t - k J dx J dyG(h j7 y - Cj)T yy (34) 



+k / dyG(hj,y-Cj)[T x (x + a,y) -T x (xi- P, y)] 



> / dyG(hj,y - Cj)u(y)[T(x + a,y) -T(xi- P,y)]. 



We dropped the integral of f(T) on the right side which resulted in the inequality in (U). The reason 
that our averagings in x and y are different is that while the width hj is a prescribed quantity, we 
have no a priori control over rjo and rji. Therefore our bounds may not involve them, and we employ 



different estimates when averaging in x. First we estimate the integral of T yy on the left side of (34) 



Lemma 2 There exists a universal constant C > such that the following estimate holds for every 
x G [x ,xi] 



« / dyG(hj,y - Cj)T yy (x,y) 



v hj 



2 C 3+ h 3 

^ J dyf(T(x,y)) + K I dy\VT\ 2 (x,y) 



(35) 



Proof. We use (|32j) to rewrite the left side of (|35|) for a fixed x E [xo,aci] as 
Cj+hj 

J dyG(hj,y- c,)T 2 



' 'T( Cj + ^) - 2T( Cj ) + T(cj - < t-Sj[T}(x) (36) 
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with 8j[T](x) = sup T(x,y) — inf T(x,y). Note that because of our choice of xq and x%, given any 



x G (xo, xi), we may find y' such that T(x, y') G [6>i, #4]. Then we may find y\, y 2 such that for any 
y G [2/1,2/2] we have T(x,y) G (#1 — £,#4 + C), and, moreover, 



Then (|ll|) implies that 



\T(x,y 2 ) - T(x,yt)\ = min(C, Sj[T](x)). 



f(T(x,y))dy > f \y 2 - yi\ 



(37) 



Applying the Cauchy-Schwartz inequality we also obtain 

^riVT^)! 2 ^ > \n*,V2)-T(x, yi )\* 
J 1 2/2 — 2/1 1 

Cj —hj 

Multiplying these two inequalities we obtain 



f(T(x,y))dy / \VT(x,y)\Zdy>^r \T(x,y 2 )-T(x, yi )\ 



\cj-hj 
,1/2. 

> ^ \T( Cj + - 2T( Cj ) + T{ Cj - hi)\ 

because of ( |36| ) and (37), and since Q < 1/2. Then ( |35[ ) follows. □ 

Furthermore, because of our choice of 770, Tji, we have for any a G [0,%] and f3 G [0,771] 



dyG(hj,y - Cj)u(y)[T(x + a,y) - T(xi - /3,y)] > (6> 3 - #2)^. 



(38) 



We use (|35|), (5|) and positivity of T 4 to rewrite (|34|) as 



C dx dy 



v hj 



K \VT\ 2 + ^f(T(x,y)) 



CA-\-hj 



+k J dyG(h j ,y-c J )[T x {x + a,y)-T x (x 1 -(3,y)]>(e 3 -8 2 )F J . 



(39) 



In order to deal with the integral term on the second line that involves T x we average ( p9|) in 
a G (0,7/o) and /3 G (0, 7/1): 



1 



rjoVi 



V0 Vl 

da J d(3- 
n 
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to get 



C J dx J dy 



T t + 



v hj 



K\VT\ 2 + ^f(T(x,y)) 

K 



c-j+h-j 



H / dyG(hj,y - Cj)[T(x + r]o,y) -T(x ,y)} 



(40) 



+— / dyG(h j ,y-c j )[T(x 1 - m ,y)-T(x 1 ,y)}>(6 3 - 6 2 )F J . 

7]l 



We bound now the term involving t]q in ( |40| ) as follows. 
Lemma 3 There exists a universal constant C > suc/i i/iaf 



- | dyG(h j ,y-c 1 )[T(x + r ]0 ,y)-T(x ,y)] <Ck(9 a -9 3 )- 2 J dy J dx\VT\ 2 . (41) 



Proof. The proof of this estimate is based on two observations. First, we have 



Cj+hj zo+r/o 



dyG(hj,y - Cj)[T(x + r/ , y) - T(x , y)] 



< 



2^ 



because of (p3|) . Second, we have 

[ dy [ dx\VT\ 2 > c(e 4 - e 3 f^. 
J J Vo 

Cj-hj xo 

This bound is established as follows. Recall that because of our choice of t]q we have 

j 

dyu(y)G(hj,y - Cj)[T(x ,y) - T(x + r] ,y)] > (0 4 - 6 3 )Fj. 



(42) 



(43) 



it 



Furthermore, recall that - — < 1^(^)1 < IMIooj on the interval (cj — hj,Cj,+hj) and thus (^3| 
implies that 



dyG(hj,y - Cj)\T(x ,y) - T(x + r] ,y)\ > 



Then we obtain using the Cauchy- Schwartz inequality: 



(0 4 - e 3 )Fj 



u 



oo,j 



c(e 4 - e 3 )h 3 < 



/ Cj+hj 

J dyG 2 (hj,y - Cj) 

\Cj-hj 

2\ 1/2 
xo+VO ~ x 

dy J dxT x (x,y) 

X{, 



V2 / Cj+hi \ 1 2 

dy\T(x ,y) - T(x + 7] ,y)\ 2 



< Ch 1 / 2 



1/9 

710 l ^ l dxT ^ X ' V ^ 

\ Cj-hj x 
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and fl42|) follows. □ 

A bound similar to (41) holds for the integral involving i]i in (|40|). We use these two estimates 
in (Mj to get 



J dx J dy 



v hj 



K\VT\ 2 + ^>f(T(x,y)) 

K 



Cj+hj Xl 



+ k J dy J dx\VT\ 



(44) 



Cj+hj 

>CFj>C' [ dy\u{y)\. 



A similar estimate holds also for the intervals Ij, on which u(y) < 0. The only difference would be 
that at the first step of obtaining the analog of (|34| ) one has to drop the integral involving T and not 
that of f(T). The rest of the estimates still hold. We use Lemma p] in ( f44|) to get after summation 



over all intervals 



-1 c 3+ h 3 



with I = k/vq. Finally, we can always add vq to the right hand side by Proposition |j. This finishes 
the proof of Theorem ||. 

3 Percolating flows 



We now consider equation 



-o 2 



T t + u(x, y)-VT = kAT + —f(T) 

K 

T{x,y,0) = T (x,y) 



(45) 



in a more general class of flows, which we call "percolating". By this we mean that there exist 
tubes of streamlines of the advecting velocity u(x, y), which connect x = — oo and x = +oo in either 
direction, as depicted on Figure |3|. We assume that the flow has zero mean ( |l^ ) and hence such 
tubes of streamlines will go in both directions. More precisely, let us assume that there exist regions 
D~j~ and DJ , j = 1, . . . N such that each of them is bounded by the streamlines of u(x, y), and the 
projection of each streamline of u(x,y), contained in either or DJ , onto the x-axis covers the 
whole real line (these projections need not be one-to-one, however). We denote by D± the union of 
all respectively. 

We will further assume that the velocity u(x, y) is periodic in space. Then it is known || |3l], |32j 
that for ignition nonlinearity ([n]) there exist periodic traveling fronts. They have the form T(x — 
ct,x,y) and are periodic in the last two variables and monotonically decreasing in the first one. 
These solutions satisfy our main condition 



dT 
~dt 



(x,y,t) > 0. 



(46) 



Our results may be generalized in a straightforward manner to non-periodic percolating flows as 
long as initial data satisfies ([l|) : 



K AT + -±f (T ) 

K 



u ■ VT > 0. 



(47) 
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Figure 3: Streamlines of u(x,y). 



However, we restrict our attention to periodic u(x,y) to simplify the presentation. 

We assume that the streamlines in are sufficiently regular, so that inside each there exists 
a one-to-one C 2 change of coordinates (x, y) — > (p, £), such that p is constant on the streamlines, while 
£ is an orthogonal coordinate for p (with a slight abuse of notation we shall use the same notation 
(p, £) in all Dj , although these coordinates may not be defined globally). Moreover, u • V£ > in 

, while u ■ V£ < in each DJ. The variable p varies in the interval [dj — hj,Cj + hj], while £ 
varies in (— oo, oo) in the set D - . See Figure ||| for a sketch of coordinates (p, £). The square of the 
length element inside each set Dj is given by 

dx 2 + dy 2 = E 2 (p, Odp 2 + E 2 (p, Odf. 



We assume that the functions E\ 2 are bounded from below: 



C^KE^p,® 



uniformly on all D^ 1 . Moreover, the function 



satisfies the following bounds: 

C- 1 <\u(p,0\<C, 



u(p,0 



E2M 



c 

< —r on Df respectively, 



(48) 



(49) 



(50) 
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Pi 

Figure 4: Curvilinear coordinates (p, £). 

with being the absolute value of the difference of the values of p on the two components of the 
boundary dD^ (recall that jD- are bounded by two streamlines of u(x,y)). Finally we assume that 
the flux density E\u does not oscillate too much on the set Dj. 



1 1 E u\\ 

— 1 — < Ei(p, d)\u(p, 0)\ < \\Eiu\\oo,j for (x,y)eDf, ||£'iti|| 00ij = sup \Ei(p,£)u(p,£)\. 

1 " (p,C)eD± 

(51) 

Note that Ei(p,£)\u(p,t;)\ is independent of £. In particular if 

u(x, y) = UHV^ = UH(—, - — 



dy dx J 

with |V^| < C/H we may choose p = H^(x,y) so that E\ = jy^W\' Then we have E\\u\ = U, so 
that ( |5l| ) holds automatically and (48) also holds for E\. Other conditions on the streamlines may 
be also easily restated in terms of the stream function ^f(x,y). 

We do not make any assumptions on the behavior of the streamlines of u(x, y) outside the regions 
D + and L>_. In particular, there may be pockets of still fluid, streamlines may be closed, etc. (see 
Figure |||) . 

Then we have the following Theorem. 

Theorem 4 Let T(x, y, t) be a solution of with the boundary conditions ([|) and fy) or 



with the initial data To(x,y) satisfying (4j_), and nonlinearity f(T) satisfying f\Q). Let each of the 



sets be of the form = {p € [cj — hj,Cj + hj]}. Then under the assumptions ft4$ ) and (|5< 
on the streamlines of the flow u(x,y) G C 1 (0), we have 

n*)>CU + E( 1 + ^) / \u(p,0\Ei(p,0^\ (52) 

for all t > 0. Here I = k/vq and the constant C in (5i) depends only on the function f(T) and the 
constants appearing in ( f^5[ ) and (5C). 



Note that the integrals on the right side of ( p2[) give fluxes of u(x, y) through the tubes of the 
streamlines. As in the shear case the pre-factor (1 + l/hj)^ 1 agrees with the homogenization limit 
ill- 
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Recall that traveling fronts for periodic flows have the form 



T(x,y,t) = U(x - ct,x,y) 



(53) 



with the function U(s,x,y) being periodic in the last two variables. It was shown in (3|, 31, 32 1 that 
in the ignition nonlinearity case (|n]) such traveling fronts exist and their speed c* is unique. Their 
existence for the KPP nonlinearities (^) was shown recently in H with c > c u , c u being the minimal 
traveling front speed. The following analog of Corollary [l] holds for percolating flows, which we 
formulate separately for the convenience of the reader. 

Corollary 2 Let f(T) be either of the ignition nonlinearity type or of the KPP type (fy). Let 

also U(x — ct,x,y) be a traveling wave-type solution of (41), periodic in the second two variables. 
Then there exists a constant C > that depends only on the function f and on the constants 
appearing in fifty and $5(\ ) such that 



( 



c>C 



1'0 



V 



Proof. Corollary § follows from Theorem |] as follows. Let U(x — ct, x, y) be a periodic traveling 
front solution of fl45| ) such that 

U(s,x + L,y) = U(s,x,y) 

and let r = L/c. Then we have 



V(t)dt 



dt / (— cU s (x — ct,x, y)) 



dxdy 
H 



dt' / [-U s (x-t',x,y)\ 



dxdy 
H 



dt' 



Q 
L 




d 

d:X 



U(x - t',x,y) -U x (x- t',x,y) 



dxdy 



H 



dxdy 
H 



c 

C+ L 



dxdy 
H 



[U(x, x + L,y) - U(x, x, y)] = c. □ 



Remark. We note that Corollary ^ implies a lower bound on the effective diffusivity in the 
homogenization regime. Recall that solutions of the advection-diffusion equation ([45]) with f(T) = 0, 

and with advection of the form u(x, y) = 

solution T of the homogenized problem 

. d 2 T 



x y 
e e 



with v(x, y) periodic, converge as e — > to the 



dT 
~dt 



13 dxidxj ' 



T(x,y,Q) =T (x,y), (xi, x 2 ) = (x, y). 



The effective diffusivity k* is a complicated functional of advection Uv(y). Explicit bounds on k* in 
terms of the magnitude U of advection are easy to obtain in the shear case, when k* may be found 
explicitly and effective diffusivity in the direction of the flow k* x ~ U 2 . Using results of [18| one 
may deduce from Corollary || that this bound applies also to periodic percolating flows, despite the 
fact that no explicit expression for k* is known in this case. 

Corollary and the stability results of [ 32 1 imply the analog of Theorem ffl. 
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Theorem 5 Let the initial data Tq(x,u) for the equation MQ ) satisfy the decay to 0,1 conditions 
1^) and fffl). Then both for KPP and ignition non-linearities we have 



lim (V) t > C 



t—>oo 



( I X -1 C 3+ k J \ 



where constant C depends only on the non-linearity f and on the constants appearing in [J^t ) and 



Now we turn to the proof of Theorem |4|. The proof is a modification of the proof of Theorem 
[|. We will again utilize Lemma [l| as well as averaging along the streamlines of u(x, y) in order to 
bound the arising averages of AT in terms of integrals of f(T) and |VT| 2 . The additional technical 
difficulties are due to the fact that two natural geometries of the problem - streamlines for the 
advective term and Euclidean coordinates for Laplacian - are in harmony in the case of shear flows, 
but at odds in the case of more general percolating flows. Moreover, while in the case of shear flows 
we gave the proof that we felt was simplest, here we will use the approach which is slightly more 
involved; however, it is better adapted for the application to cellular flows in the following sections. 

Let us consider a region DJ = {(/>, £) : p € [cj — hj, Cj + hj]} with u ■ V£ > 0. Introduce notation 
kj(p) = G(hj,p — Cj)Ei(p, £)| (kj does not depend on £ by incompressibility of u(x,y)) and, 

similarly to the shear case, 

c i+ h i 

Fj= J k J (p)dp. 

We can find £i > £o such that 

J k j (j>)T(p,£ )dp = 9 i F j , J k J (p)T(p,Ci)dp = e 1 F J , 

Cj—hj cj—hj 

and for every £ G [£o>£i] we have 

Cj+hj 

0iFj< J kj(p)T{p,(i)d P <e±Fj, 

Let us denote by DJ the region bounded by the curves p = Cj ± hj and £ = £o,£i- This region 
depends on time, but we will suppress this dependence in notation. Theorem ||| will follow from the 
following 

Theorem 6 Let T(x, y, t) be a solution of ([^j with the boundary conditions (^) and fy) or 
with the initial data To(x,y) satisfying ftjlj), and nonlinearity f(T) satisfying ^). Then under the 
assumptions and $5(\) on the streamlines of the flow u(x,y) € C 1 (r2), we have for every time t 

T t + k\VT\ 2 + ^/(T) ) dxdy > C((, /„) (l + ~) ' (0 4 - #i) 3 ^ 




with the constants fo and #^4 defined by (11) 
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Remark. A similar lower bound on the integral over Dj is easier to obtain, and also suffices to prove 
Theorem ^. We chose, however, to formulate Theorem ^ in this stronger version since this is what 
we will need when dealing with cellular flows. 

Proof. Let #3 = #4 — (64 — #i)/3, and 9 2 = 61 + (#4 — #i)/3. Then we may choose 770 and 771 similarly 
to the shear case, namely 

! Cj+hj 
r ? >0: J dpk j (p)T(p,£ + r))=9 3 F j 



m 



1 1 j 

inf|r ? >0: J dpkj(p)T(p, £1 - 77) = 9 2 Fj 



We choose a G (0, 770) and f3 G (0, 771), integrate (^) over £ G (770 + a, r\\ — f3) and average in p with 
the kernel G to get 

/ di J dpG{hj,p- C j)E x E 2 T t - K J d( J dpG{hj,p- C j)E x E 2 [T xx + T yy \ (54) 

Cj+hj 

> [ kj^mp^o+^-Tip^t-pjjdp 



> 



F;. 



In (|54|), we dropped the term involving /(T) on the right-hand side, as we did in the shear case. We 
now look at the term involving the Laplacian: 

j d£ j dpG{h j ,p-c j )E 1 {p,0E 2 (p,0[T xx {p,0+T yy {p,0} 

§0+Q c j—hj 

= Y.j d5 J J d P E 1 (p,OE 2 (p,0[T xx (P,0+Ty y (p,0]- 
J 5o+o Cj— S 

Let us denote D$ a p the region bounded by coordinate curves p = Cj ±5 and £ = £0 + 01, £ = £1 — /?. 
Using Green's formula, we rewrite the last two integrations for fixed 8 as 

6-/3 Cj +<5 

/ / dp J Bi(p,^2(p,0[^(p,6+^(^C)] = 



/ / ATdxdy 



dD Sa « dn 



dT , 
-—ds 



a-/3 



<9T dT 

uj(cj + 5^) — {cj + 8,0 -uj(cj -8,Q — [cj -8,0 



dp dp 

r r <9T or 

+ y dp uj-^p,^- (3) — (p,^- p)-uj- 1 {p,^ + a)-—(p,^ + a) 



(55) 
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We used the definition ( [49|) of the function uj(p,0 = E 2 (p,0 1 Ei(p,0 in the last step. The average 
of the term on the first line in (|55|) is estimated by the following Lemma, an analog of Lemma 



Lemma 4 There exists a constant C > such that we have for all £ G [£o,£i] 

hi 



K 
kj 



dS 



dT dT 
uj(cj + ^0-r^( c j + S ,Q ~ w ( c j - ^0~Q^( c j ~ 8,0 



v hj 



Cj+h 



d P E x E 2 f(T)+K / dpE x E 2 \VT 



where the constant C depends only on constants in the bounds (5C) and (48). 



Proof. We will show that 

f dT 
J ddu( Cj + 5,0-g-(cj + 5,0 



< 



cf 1/2 C-> 



Cj+hj 



Cj+hj 



(56) 



k J d pT 2 p (p,0 + ^ J dpf(T(p,0) 



(57) 



for all £ G [£o, £i]- A similar estimate holds for the second term on the left side of (p6[). By definition 
of £o an d £i for every £ G [£o,£i] there exists po G (cj — hj,cj + hj) such that T(p ,0 G [^1,^4]. 
Then given £ G [£o>£i] we have two possibilities. First, assume that T G (6\ — £,#4 + C) for all 
p G (cj — hj , Cj + hj). Then we have for such £ 



I dpf(T(p,0)>Cf J dpuj 2 (p,0, 



which implies that 

h, 



c-j+hj 



1/2 



J dSu;(c J + 5,0-^(^ + 6,0 <\ J dpu\p,0 

( c ]+ h j 



I Cj+hj 



1/2 



d P T 2 Ap,0 



<Cfo 



-1/2 



1/2 /cj+hj 



1/2 



V 



dpf(T(p,0) 



d P T 2 Jp,0 



< 



Cf, 



-1/2 



^3 

i+hj 



k / d pT 2 (p,0 + ^ / dpf(T(p,0) 



The other case for a given £ is that the temperature T(p,0 drops out of the range (#1 — £,#4 + 
for some p. Then we may find p' G (cj — /ij, Cj + /ij) such that T(p', = #1 — C or ^""(/A £) = $4 + C> 
and, moreover, T G (#1 — £, #4 + C) f° r an P between p' and po- Then we have \T(p' , £) — T(po, £) | > £ 
and hence 

\ 1/2 /^+^- \V 2 1/2 

/ dp/(r(p,£))j I | d P r p 2 (p,£)j > ^^X-^/oip'-poij > C /V 2 . 
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Then the estimate ( |57|) also holds in that case since 



r dT 



< C 



for all £, as can be seen from integrating by parts in p and using (|50|). Since C < 1/2, (p7|) holds in 
both cases. □ 

In order to bring the term in the second line of ( |55| ) into a form convenient for analysis we 
average (|55| ) in a and 0. Let us consider an estimate on the second summand after averaging in a 
over [0, 7/0 ]; the other summand is treated similarly. 



Lemma 5 There exists a constant C, depending only on constants in bounds ( (8 ) and ( 5C ), such 
that 



Cj+8 



■ J d5 J da J dpu 1 (p, £ + a)^(p, £ + a) 



% h. 



< C 



fo^ie^-O^mm^-e^Xr 1 ^- [ dxdy 



(58) 



voiij j 



^f(T(x,y)) + K\VT(x,y)\ 



-(O^-e^K I dxdy\VT(x,y)\ 2 



D1 



Proof. We have to estimate the following expression 

hj }j c j+<5 



- j d5 j da J dpaj- l (p,£ + a)—(p,£ + a) 

cj-6 
^ h j c j+ s 

= — dS / dplu' 1 (p,£ + Vo)T(p, £0 + 770) - w _1 (p,^ )r(p, Co)] 
Tin J J 



(59) 



Cj+5 



Iff f duj^ 



The first term in ( |59|) is evidently bounded by 



1 

m 



j d5 J dp[uj 1 (p,Co + 7/o)T(p,Co + T] ) ~v 1 (p,(o)T(p,£ )} 

Cj -5 
Furthermore, we claim that 



< 



cjl 

Vo 



(60) 



dxdy\VT(x,y)\ z > 



2 C(6*4 — 3 ) 2 hj 



T)0 



(61) 
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This is shown exactly as the estimate (|4^) in the proof of Lemma || for the shear case, given the 
assumption (pC|) on uj. We combine (|60|) and (|6l|) to obtain 



1 



J d5 J dp[u 1 (p,^ + Vo) T (p,€o + Vo) ~ u 1 (p,£o)T(p,£ )] 

Cj -S 

< C(0 4 - 6z)- 2 h 3 J dxdy\VT(x,y)\ 2 . 



(62) 



« 1 



Next we estimate the second term in (|59|) : 



l 

m 



dS 

go Cj-8 



BbJ~ l 

< I dp-—(p,C)T(p,0 



< Chi 



(63) 



because of ( |50| ) . Notice that if r]Q < hj , then Lemma |5] follows directly from (^) . Hence we can 
assume rjo > hj . The following final Lemma allows us to finish the proof of Lemma [|. 

Lemma 6 Assume that r/o > hj. Then there exists a constant C, which depends only on the constants 
in bounds (5C) and (48), such that 

J f(T)dxdy J |Vr| 2 ^>C/o(04-ei) 2 (min((e 4 -^i),C/2)) 2 ^ 2 . 



Remark. This lemma is much easier to prove with Dj instead of Dj~, basically the argument of 
Proposition [l] applies. The simpler version is also sufficient for the proof of Theorem ||, but we need 
this stronger version for the proof of Theorem |(| 

Proof. Let us denote \S\ the Lebesgue measure of a measurable set S. We know that 

Cj+hj 

J kj(p)(T(p,( Q ) - T(p,^))dp = (0 4 - e 1 )Fj. 

Consider two sets 



Then 



and hence 



S x = {p\ T(p, £o) - T(p, Ci) < (0 4 - 9x)/2}; S 2 = [ Cj - hj,c 3 + hj] \ S x . 
kj(p)(T(p,£ Q ) -T(p,^))dp < (0 4 - O^Fj/2, 

k j { P ){T{p^)-T{p^ 1 ))dp>{e A -e 1 )Fj/2. 



Si 



S-2 



Since T(p, £o) — T(p,£i) < l, and kj(p)/Fj < 2h- 1 by the properties of G and ( |5l| ) we must have 
\S 2 \ > (#4 — d\)hj/A. Choose po £ S 2 such that 



6 



<%f(T(p ,Z)) < 



\S 2 



f(T)dxdy 



D J 
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and 

So D+ D + 

J 3 

with C depending only on the constants in (j5C|). Since po G S2, we have T(po,£o) — T(po 5 £i) > 
(6*4 — 6i)/2. Assume first that there exists £ 6 [£o)£i] where T(po,£) € [#1 — C/2, 64 + £/2]. Then we 
can find £ < 6 < 6 < 6 such that for every £ G [£ 2) 6], T (Po,€) G [#i - C, #4 + C]> and 

r(p ,6) -T(po,&) > min((0 4 - 0i)/2,C/2). 

Therefore, 



6 

/ /(r)de>/o(6-6) 

So 



while 



(min((04-^i)/2,C/2)) 2 



6) 



e 3 -6 

Combining these estimates and the definition of /0o> we find 

| f(T)dxdy J \VT\ 2 dxdy > C/ o (min((0 4 - #i)/2, C/2)) 2 1^ 2 , 



which proves the lemma in this case. 

The other case we have to consider is that for every £ G [£o>£i]> T(po,£) £ [6\ — C/2, #4 + C/2]- 
Then we will be able to find a point £ G [£o>4"i] so that there is a drop in temperature along the 
curve £ = const. Assume T(po,£) > #4 + C/2 for all £, the other case being similar. By mean value 
theorem, we also have that for every £ G [£o>£i] there exists p\ such that T(pi,£) G [#1,6*4], and so 
T(pi,£) < 6*4. Similarly to the above, we can find £ such that 

Cj+hj 

I dpf(T(p,i))<-^— [ f(T)dxdy<- [ f(T)dxdy 

J si so J m J 

Cj-hj I)- i) 

and 

Cj+hj 

[ dpT*(p,i) < [ T 2 p dpdi < - [ T 2 dpd£ < — [ \VT\ 2 dxdy. 

J p £1 - £0 J y Vo J Vo J 

0,-hj d+ d+ 

An argument identical to the one we used in the previous case establishes that 

Cj+hj 

f(T)dp > f ( Pl - p 2 ), 



while 

Cj+hj 



Pi) 
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for some Cj + hj > p\ > p2 > Cj — hj . This implies that 

J f(T)dxdy J \VT\ 2 dxdy>CM 2 V 2 >Cf ( 2 h 2 . 



5 / 



Thus Lemma H is proven. □ 

Lemma || and the assumptions (HSI) and (|50|) imply that 



dxdy 



^/(D + ^IVTI 



> C/ 1/2 (^4 - 0i)min((0 4 - ^i)/2,C/2)^. 



Therefore we deduce from ( J63D that 

N £o+»?o c J+< 5 



1 



Vohj 



BbJ~ l 

d5 j d£ j dp—(p,QT(p,t) 

fo 



(64) 



< C/ 1/2 (^ 4 - flirWtt^ - 0i)/2, / 

D+ 



dxdy 



-f(T) + k|VT| 



Then (^) , (^) and (|64]) imply Lemma |5[ □ 

Averaging of (pq) in (5 G [0, &i] produces terms similar to the left side of fl58|) , which are bounded 
by the same quantity. Now we put together our estimates on Dj , that is, equations (j54l) and (55), 
and Lemmas |I| and [| : 



y dxdyG(hj,p-Cj)T t + CC 1/2 ^~ ^ 



Dl 



Dl 



C/o^^-SirW^-e^C)- 1 ^ J dxdy 

3 D+ 

j 

+C{9 4 -e 1 )- 2 K J dxdy\VT\ 2 >C(0 4 -9i)Fj. 



—f(T) + k\VT\ 



^f(T) + K \VT\ 2 

K 



+ 



Dl 



Therefore we obtain 



1 + 



vohj 



dxdy 



D 7 



T t + —f(T)+K\VT\ 2 



>C(C,/ O )(04-ei) 3 ^ 



3 3 

>c(c,/ o )(04-0i) 3 J e^ohp^m, 



which proves Theorem ^. □ 

The estimate on the regions DJ with u ■ V£ < is similar. The only essential difference is that in 
inequality (|5^) we drop the term involving Tt but keep the one with f(T), which does not make any 
difference in the final result. Then Theorem [| follows after summation over all Df fr om Lemma Rl 
and Proposition |]. □ 
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4 Cellular flows: the main result 

Now we consider (|2|) in a cellular flow: 

T t + u- VT = kAT + ^f(T), (65) 

K 

that is, flow with closed streamlines, and establish a lower bound for the burning rate. For simplicity, 
we limit our consideration to one typical representative class of cellular flows, given by the stream 
function 

ip(x, y) = UH sin — sin (66) 
H H 

on the strip (— oo, oo) x [0, ttH] (for convenience from now on H will be the width of the strip divided 
by 7r). The flow u(x,y) is given by 

u(x, y) = UV ± ip(x, y) = U ( sin — cos — , — cos sin — ) . (67) 

The streamlines inside a period cell are depicted on Figure pi. The results we prove can be extended 
in a direct way to the periodic cellular flows of more general form. 
We will further assume that the Peclet number is larger than one: 

Pe = — > 1 (68) 

K 

and the size of the cell is larger than the laminar front width: 

H^ 1 = V„ (W > 



Conditions (68) and ( p9|) are natural for flows of large amplitude and for thin fronts. Moreover, our 
results may be easily adapted to the other regimes, where (38) and ( |69[ ) are violated. We introduce 
also the turnover time t u , and the chemical reaction time r c : 

H K 

T u = — , T c = -g. (70) 

It turns out that the ratio t u /t c is the crucial parameter for burning in the cellular flows. 
Finally we assume that T(x, y, t) satisfies the usual boundary conditions and that 

dT 

-^(x,y,t)>0. 

As we noted previously, this condition is satisfied as long as it holds initially. 

Theorem 7 Let T(x,y,t) be a solution of (21) with the boundary conditions (W and either (Qj or 
(£|), and the cellular flow given by (j^P- Let the initial data T (x,y) satisfy $24), (^) and (Q), and 
let the non-linearity f(T) be of either ignition or general KPP type. Furthermore, assume that (|#£p 
and M3) hold. Then we have for any time t 



[Cu ^ + C 2 )v , if t c <t u 
V(t) > \ f f x 1/5 \ (71) 
1 Cl(^) +C 2 J^o, if t c >t u 

The constants in the inequalities depend only on the reaction f , more particularly on constants /o, 
£, and 9^ — 0\ that appear in ^L 
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Remark. As we noted above, in order to avoid excessive details, we chose not to formulate Theorem 
in the exhaustive form which goes through all possible relationships between parameters (large k 
limit, small H limit, small vq limit). The reader will find it not difficult to extend the results we 
prove to the above mentioned regimes. Theorem ^ is formulated here for the range of parameters 
that appears to be physically reasonable for most problems of interest. 
Furthermore, we have the following corollary. 

Corollary 3 Let f(T) be of ignition nonlinearity type /f7Z[), or of the KPP type ^) and let c be 

the speed of a traveling wave-type solution T(x,y,t) = U(x — ct,x,y) of ffih]), periodic in the second 
two variables. Then there exist constants C\ t 2 > which depend only on the function f and on the 
constants appearing in ( f^ ) and (5C ) such that 



c > < 



CiJ^ + C 2 )v , if r c < t u 

1/5 



Corollary || follows from Theorem since the traveling front profile U(s, x, y) is monotonically 
decreasing in s (see || for the ignition case, and Q for KPP). Then Theorem || follows immediately 
from Corollary ||, results of |3l| and the argument we gave in the proof of Theorem [l]. 

The proof of Theorem [j] is a boundary layer argument that proceeds, roughly, as follows. The 
temperature drops from one on the left to zero on the right. We will watch the temperature in the 
layers of width h formed by streamlines near the boundary of the cells. The drop of temperature 
in these layers may occur inside the cells or over the diffusive interfaces. The first estimate, which 
we call advective, shows how much the cell must contribute to the bulk burning rate if a certain 
drop of the temperature (in the range [#i,#4]) takes place along the streamlines inside the cell. It 
is reasonable to expect that the drop over the cell will be small when advection is strong since it 
mixes the fluid inside the cell quickly; in an analytic form this intuition will translate into a large 
lower bound for the burning rate if the temperature drop is significant. The second estimate, which 
we call diffusive, gives a lower bound for the burning rate given certain drop of the temperature 
between the two cells. We do expect the temperature to drop on the boundaries, and hence the 
lower bound is only effective if we choose h in an appropriate way, sufficiently small. Finally, we 
prove the reaction estimate, which takes into account the total area of the region over which the 
temperature drops. These estimates will be brought together to establish the lower bound for the 
bulk burning rate using an appropriate optimization argument. 



5 Cellular flows: regularity of the streamlines 

Our first objective is to define appropriate curvilinear coordinates on the cells, and to show that 
these coordinates satisfy certain technical assumptions that we will need. The natural choice of the 
coordinate p, which is constant on the streamlines, is 

ip{x,y) . x . y 

P(x,y) = — - — =Hsm — sin- . (72) 

We have certain freedom in the definition of the orthogonal coordinate £ along the streamlines : 

V£ = QV x p, (73) 
where Q is some function which should satisfy 

Vp • S7Q = -QAp. (74) 
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It is easy to compute that 



dp = x p d x + y p d y = -^j-d x - ^d. 



where J = p x £ y — p y £ x = — Q|Vp| 2 . Hence 



0, 



1 



Vp-V, 



J 



and so we have from (f74|) 



P |Vp| 2 
dQ QAp 



dp \v P \ 2 ' (75) 

We will choose Q so that Q(H/2,£) = 1; equation ( |75| ) then allows us to define Q in the region 
H > p > : 



Ap(fe,Q 



/ |V P (h,5)| 2 



(76) 



We have the following auxiliary 

Lemma 7 In the region H/2 > p > we have the following bound for the function Q : 

< e" 1 < Q{p,0 < e. 
Proof. Consider formula (|76|). We have 



and 



|Vp| 2 



= -j^p(x,y), 



x . y\ 



x y 

cos — sin ) +| sin — cos — 



— sin — 
H H 



V ^2 



2P 2 



in the region p < H/2. Hence in the region H/2 > p > 0, 

|Ap| 



(77) 



|Vp|= 



< 2/ff. 



Therefore, by (fTq), we have e 1 < Q < e in this region. □ 
Recall our notation 

(fa; 2 + cfy 2 = Efdp 2 + E%df, u) = ^. 

The next proposition summarizes some of the properties of the coordinates (p, £) in a region of 
interest to us. 



Proposition 2 For the cellular flow defined by ^6% ) and coordinates p, £ defined by (7i), fj7jj), and 
(7t), the following bounds hold in the region H/2 > p > : 

Ei, 2 (P,0>C (78) 

< C- 1 > uj(p,0 < C (79) 

^(p,^<CH- 1 (80) 

8oj 



(p,0 



^cir^iogOo/fl-)!. 
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Proof. Direct computation gives that E\ = l/|Vp|, Ei = l/Q\Vp\. Then Lemma ^ and the fact 
that |Vp| < 1 imply (|78|). It follows from the definition of uj that to = l/Q, and hence Lemma [?] 
implies (M). Next, 



proving 



Notice that 



Finally, 



duj 



0( 



dp 
2 



1 [Agj 2e 
Q|V/9| 2 - if 



Qif 



/ 1 



H/2 H|Vp| 2 



9„. 



y QWp? x Q\v P \ 2 V 

A straightforward computation using © and (0) leads to 



<9u 



< 



c_ 

H 



p dh 

H/2 h 



<CH- l \log(p/H)\. □ 



6 Cellular flows: Advective estimate 

Let us introduce some notation. Within the cell, we will normalize £ by letting it be zero in the 
negative direction of the x axis (assuming that the origin has been placed in the center of the cell) . 
We will denote by L the value of £ in the positive direction of the x axis. In every cell, we will 
consider a tube of streamlines bounded by p = h and p = 3h; h will be always assumed to be less 
than H/6. We set 

k(p) = G(h,p-2h)E 1 (p,0H P ,0\- 

The fact that k(p) does not depend on £ is a direct corollary of incompressibility of the flow. 
Moreover, with our definition of p for the cellular flow, we have k{p) = UG(h,p — 2h). We also 
denote 



3/i 



F 



k{p) dp. 



As a corollary of Proposition ||, in the strip h < p < 3h we have all conditions on E\^ and uj that 
have been necessary for the percolating flow estimates. In particular, we have 



< CH- 1 |log(p/F)| < ChT 1 , 



(81) 



since h < p < H . We first state an estimate very similar to Theorem ||. 



Theorem 8 Assume that within one cell C, there exist two values £0, £1 such that for some sq, s\ € 
[01, #4], sq > s\, we have 



3h 



3h 



J k(p)T(p, Co) dp = s F, J k(p)T(p, £1) dp = siF, 
h h 

3h 

siF < J k(p)T(p,Odp < s F for £ G [f ,£i]. 
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Let D be the region bounded by the curves p = h, p = 3h, £ = £o and £ = £i- Then 

T t + k\VT\ 2 + ^f(T)) dxdy > C(C, /o) (l + £) ' (*> - ^i) 3 ^ (82) 



The constant C in (lei) depends only on parameters £ and /o 0/ i/ie reaction f , that appear in (11 ), 
and on the constants in the bounds of Proposition |j. 

Proof. The proof is exactly the same as for Theorem || for percolating flows. We need only to 
replace 6\ with s\, 64, with so, and set Cj = 2h, hj = h. □ 

The estimate ( |82[ ) works well if there is a significant change of the temperature along the stream- 
lines within the cell. But it has a serious flaw if the temperature drops gradually and there is little 
change of temperature inside any cell. The factor (sq — s\) 3 on the right hand side of the estimate 
makes it rather inefficient. Our next goal is to derive an estimate which has linear dependence on the 
temperature drop when the drop is small. Our main measurement tool for the temperature within 
the cell will be the following average 

Co+A 3h 

{T) ^ = 2lF J di j k (p) T (P^)dp- 
to— A h 

We will normally take A so that the region of averaging in £ covers about half of the width of the 
strip when £0 = or £0 = L, hence A ~ H. Our measure of the temperature change along the cell 
will be the difference between such averages for different £o- 

Theorem 9 Assume that for all £ in a given cell we have 

3h 



h 



Then for any £0, £1 we have for A> h 

£i+A 3h 



J J (r t + ^|VT| 2 + ~f(T) \ E 1 E 2 d£dp > C(/ , C)(l + ^ 

A I, \ / 



\(T) i0 -(T)^\F, (83) 



Zo-A h 

where C depends only on /o, £ and constants in the bounds of Proposition || 

Proof. We will consider the case where {T)^ > (T)^, the other case being similar. Let us integrate 

v 2 

T t + u-VT- kAT = —f(T) 

in x and y over the region where £ varies from £2 to £3, where £2 G [£o — A, £o + A], £3 S [£1 — A, £1 +A], 
and p takes values between h and 3h, with kernel G(h,p — 2h). We obtain 

£3 3h £3 3h 

J d£ J d P G(h, p - 2h)E l E 2 T t (p,£) - k J d£ J dpG{h, p - 2h)E 1 E 2 AT 

?2 h t% h 

3h 

> I dpk(p)(T(p,&)-T(p,Z 3 )). (84) 



h 



30 



As usual, our goal is to control the Laplacian term. Rewrite it in a familiar form 




dT dT 
—u(£,2h + 5)-—u{£,2h-d) 



2h+S 



+ 



dp 



2h-8 



(85) 



We carry out two more averagings in (]8 



AA 2 



<%2 



$ -A fr-A 

to be able to estimate the second part of the Laplacian term (|85|). Let us start by estimating the 
first expression in square brackets, more precisely, 

h ^ 



o 6 



(the other part is estimated similarly). Notice that by the assumption of the theorem for every £ 
there exists p (depending on £) such that T(p,£) € [61,04]. The estimate of this term now follows 
step by step the estimate of the same term in the proof of Theorem |(| We get that for every 



h 3h 



K |VT| 2 + —f(T) 

K 



E x E 2 dp. 



(86) 



The integrations in £2,3 simply average out, they are not needed for this term. 

Let us now consider the estimate of the second expression in square brackets in fl85"D, more 
particularly the first summand (the second one is estimated in the same way). Averaging in £2 
simply disappears since there is no dependence on this variable, and we are left with 



h 2h+6 £i+A 



2hA 



Jd5 J dp J d^^i^p). 



2hA 



1/2,-1 



hv 



2/i-5 

The following lemma is the crucial step in the proof. 
Lemma 8 Under conditions of Theorem |^ we have 

h 2h+8 £i+A 

I d5 

2h-S £i-A 

where the constant C depends only on the bounds in Proposition |^. 
Proof. Integrating by parts and using (|8l|), we find that 

Si+A 



dT 

dp I ^3— cj _1 (^3,/j) 



D 



K |VT| 2 + —f(T) 

K 



dxdy, (87) 



r dT 

j dts—u-H^p) 



31 



and therefore 



2hA 



h 2h+8 £i+A 



h 



<Ck[- + 1) <Ck. 



(89) 



Jds J dp J d^—u-H^p) 

2h-S %!-A 

We have to consider several cases. 

Option 1. There exists po such that for all £ E [£i — j4,£i + A], we have T(po>£) ^ [$i — C/2, #4 + C/2]- 
By mean value theorem, for every £, there also exists /?i(£) such that T(pi,£) E [#i,#4]. Then 
for every £ we can find /?3 > ^2 such that |T(p2>£) — ^G°3>£)l = C/2, and for every /? E [/02,P3], 
T(p, £) E [#1 — C/2, #4 + C/2]- Therefore, as we have seen before, for every £ E [£1 — .A, £1 + A] we 
have 

3fc \ V2 , 3/l , 1/2 



dpT 



dpf(T) 



> 



C/o 



1/2 



Hence, using Cauchy-Schwartz and integrating in £, we get 

37i (i+A 



Oil ilT" / 2 \ 

J dp J dmE 2 L\VT\ 2 + ^f(T))> 



h 



vo(fo /2 A 



From this inequality and ( p9|) our lemma follows since A > h. 

Option 2. For every p, there exists £4 such that T(p, £4) E [#i — C/2, 6*4 + C/2]- Here we have to 
consider two distinct sets of p. First, assume that T(p, £) E [#i — £,#4 + C] f° r all £• Denote the set 
of all such p by 5i . For p E Si we have 



/ 

€1-^4 



| /(T)d£>C/ | ^ 2 (p,£)d£, 



and so 



6 +-4 

/ 

6 -A 



1/2 



<Cf, 



-1/2 



/(T)d£ 



6+A 

/ 

\Hi-A 



IT" ' 

/ 7fd£ 



1/2 



(90) 



< cf^W 



K |VT| 2 + ^/(T) 

ft 



£i£ 2 d£- 



The second case is that there exists £4 E [£1 — A, £1 + A] such that T( / o, £4) ^ [0\ — £, #4 + C]- Denote 
the set of all such p by S%. For p E S2, we can find £5 and £6 such that \T(p, £5) — T(p, £g)| = C/2, and 
T(pX) G [#i — C>#4 + C] f° r all £ between £5 and £g. In this case, similarly to the above reasoning, 
we have 

Zi+A £i+A 
J f(T(p,0)d£ J T|(p,£)d£>CC 2 /o, 



6 -A 



and so 



£i+A r 



/ - 1/a rV 



kIVTI 2 + —f{T) 



6 -A 



^1^2 <i£ > C. 
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Figure 5: The regions D n . 



Then (gg) implies that for p € S2 we have 



dT 

^3-^-W _1 (^ 3 ,p) 



K\VT\ 2 + ^f(T) 

K 



E^dt. (91) 



Combining the two estimates ( p0| ) and (91), and taking into account that A > h, we obtain the 
result of Lemma ||. □ 

Now we can finish the proof of Theorem [|. Taking into account ( |86|) and (|87|) we see that 



3/1 



5o-A ft 



J d£ J dpG(h,p-2h)E 1 E 2 T t (p,O + Cf 1/2 (- 1 : ^ J dtJd P E 1 E 2 

-Ah 

to+A £i+A 3/1 

j2 f J d€ 3 Jdpk(p)(T(p,t 2 )-T(p,{ 3 )), 



K |VT| 2 + ^/(T) 

K 



> 



AA 2 

io-A £i—A h 
and this implies 

J J '(T t + K|VTj 2 + ^f(T)\E 1 E 2 dtdp> C(/o,C)(l + ^) |<T> & -(T) & |F. □ 

7 Cellular flows: Diffusive estimate 

Our goal in this section is to estimate the burning rate from below in terms of the jump of the 
temperature across the interface separating two cells. We will only consider this estimate in the 
context of the particular cellular flow we are studying, though it can be easily extended to a more 
general situation. Consider two neighboring cells, which we denote C n -\ and C n . We will look at 
two regions D 2n -\ C C n _i and D 2n C C n which are symmetric under reflection with respect to the 
line separating the cells (see Figure ||). The region D 2n is bounded by the curves p = h, p = 3h, 
£ = —A, and £ = A. For simplicity, we choose p and £ coordinates in cell C n -\ so that p changes 
from —3h to —h and £ from —A to A on D 2n _\ (see the figure). Notice that this is different from 
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the choice of coordinates we employed locally in each cell in the previous section; in terms of the old 
coordinates, T>i n -\ C C n -i is bounded by the curves p = h, p = 3h, £ = L — A, £ = L + A. With the 
new choice of coordinates, we denote by I n the part of C n _i U C n bounded by the curves p = —3h, 
p = 3h, £ = —A, and £ = A. We also denote, in the extension of notation of the previous section, 



(T) 



1 



2n-l 



2AF 



T(p,Z)k(p)d£dp, 



D 2 „- 



(T) 



2n 



2AF 



T(p,Z)k(p)dZdp 



(92) 
(93) 



D 2 n 



(k(p) is defined in a natural way as k{—p) when p < 0). 



Theorem 10 There exist constants depending only on the constants in the bounds of Proposi- 
tion such that 

k J | VT| 2 dxdy > ((r)^ - (T) 2rt _ 1 ) 2 (94) 



for A = C 2 H. 

Proof. Notice that 




-3h 



k(p)T(p,i)d P 



k(p)T(p,0dp\ < 



A 3h 



A~3h 



1/2 



dT 
dp 



(p,0 



dp 



<Cj^\j \VT\"dxdy 



3h 



3h 



2AF 



dH / k(p)dp 



h 

1/2 



-3ft 



dT 
dp 



(7,6 



Comparing the quantity we estimated with the definitions (92) and (pg), we see that ( |94|) follows 



□ 



8 Cellular flows: Reactive estimate 

It is worthy to note that the diffusive estimate we proved in the last section is quadratic in the drop 
of the temperature over the interface. With the estimates we currently have we could not prove 
any lower bound for the burning rate, since the temperature could stay constant inside the cells (so 
that advective estimate does not give us anything) and drop in extremely small increments over the 
diffusive interfaces. There is no lower bound for such scenario because of the quadratic dependence. 
But in this case, the region where T S [#i,#4] would be very large, and hence we can hope to get 
a lower bound on / f(T). The following theorem is a rigorous expression of the above idea. We are 
mostly interested in the case H > k/vq but include the case H < ^ for completeness. 

Theorem 11 Assume that in a given cell C n , 

3h 

1 jT(p,t)d P z[e x M 
h 
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for every £. Then 



n 



k\VT\ 2 + ^ 



f(T) 



dxdy > C(fo,C)voHmm(l, 



(95) 



Proof. The proof uses an argument we already used several times when proving advective estimate. 
Consider a region in C n bounded by the curves p = h and p = H/2. By mean value theorem, for 
every £, there exists po € [h, 3h] where T(po,£) £ [#i,04]. All £ fall into two cases. In the first case, 
T(p, £) e [0i — £, 04 + C] f° r every p in the region p £ [/i, JET/2]. Then 



i//2 



f(T(p,0)dp> 



hHvl 



In the second case, there also exists pi G [/i, fT/2] such that T(pi,£) £ [6\ — C,#4 + C]- I n this case 
we can find p2, P3 such that \T(p2,£) ~~ T(ps,£)\ > (, and for every p between p2 and p3, we have 
T(p,£) € [#i — Ci ^4 + C]- I n this case, by the usual argument, we have 



/H/2 



fo /2 C < 



H/2 



1/2 



H/2 



f(T)dp 

\h h 



\VT\ 2 dp 



< 



2v 



K |VT| 2 + %(T) 



dp. 



It is easy to show following Lemma and Proposition ^ that the length of the interval of integration 
in £, 2L, satisfies 2L > CH, where C is some universal constant. Integrating over £, we then get the 
result of the theorem. □ 



9 Cellular flows: Putting it together 

Now we have all necessary estimates to establish the lower bound on burning rate in the case of 
cellular flows. We begin with an auxiliary computation. Its goal is to choose the right value of 
h depending on other physical parameters fixed in the problem. The argument below may not be 
completely rigorous; we will fill in the gaps in the actual proof. Consider the sequence of averages 
djD), (|93|) introduced in Section [?] devoted to the diffusive estimate: 

••• , (T) 2n , (T)2n+l, (r) 2 ( n +l)> {T)2n+3, ■ ■ ■ 

where n varies from — oo to oo. The values of the averages tend to 1 as n — > — oo and to as 
n — > oo. Assume that (for sufficiently large U) the change is gradual, and there exists a number of 
consecutive numbers n where all averages lie in the interval [0i,#4]. Moreover, assume that in all 
cells corresponding to these values of n, we have 

1 [ 3h k(p)T(p,Odpe[e 1 ,e i ] 

r Jh 

for all £, so that the reactive and linear advective estimates may be applied in these cells. It does not 
concern us here that h has to be chosen yet; as we mentioned above, this is an auxiliary computation 
and the rigorous argument will appear in the proof. Let us denote 

<57n,a = ICO27H-I — (T)2n\ 
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the change of the averages controlled by the advective estimate, and 

ST ntd = \{T) 2 ( n +l) - {T)2n+l\ 

the drop of the temperature controlled by diffusive estimate. We also set 5T n = 5T nA + 8T n ^ which 
we call the total drop over nth cell. Then we have, according to the advective estimate (83), 



/ (t, + ,\VTf + £ /(T)) > cob, (1 + 



C n 



u v # 



(96) 



where the constant C depends only on /o, C an d fixed constants associated with the geometry of 
streamlines of the flow. According to the diffusive estimate (p4|), we also have 



J (k|VT| 2 + |/(T)) 



(97) 



CnUC„-|-l 

Also for every C n from the region we consider, we have by the reactive estimate (pi) that 



Cn 



VT[ 2 + ^/(T) 



^ > C(/ , C)^omin(l, ^°) = C(/ , C)«d 



(98) 



under the assumption Hvq/k > 1 (see (|69|)). Fix some cell C n , and introduce a parameter < a n < 1 
such that 

a n ST n = 5T Utd , (1 - a n )5T n = 5T n , a . 
According to (|96|), (|97|), and (p3), we have 



(1 + ^ r V aT »< 1 - < " J + T 15T » + "» 



(99) 



Over our sequence of cells, the value of temperature averages falls by a fixed amount, 64 — 61. We 
are going to assume that this fall is gradual, and consider the lower bound on the contribution of 
the cells C n and C n +i to the total burning rate normalized by the temperature falloff 5T n over the 
cell C n . Namely, denote 



CnUCn + l 

then the burning rate V satisfies 



v n = J (r t + k\vt\ 2 + ^f(T)) 



dxdy 

~1T' 



Equation (p9|) may be rewritten as 



V n > C 



1 + 



k \- l Uh 



naz 



Moreover, 



voh J H 

5>r n >0 4 -0i, 



5T n . 



(100) 
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where summation is taken over the set of cells where all estimates apply. In the lower bound fllOOD 
a n and \8T n \ are out of our control, but we may choose h optimally. Hence if we denote 

B(k,v ,U,H)= max { min ( (1 + JLyiTLn - a n ) + ^5T n + ) 1 , 

V ' ' ' ' 0</KJf/6 \o<a„<l,0<<5T n <2 ^ V h' ' h 5T n J J 

then 

V > CB(k,v ,U,H)(6 a -Ox). 
Therefore our goal is to find B(n, vq, U, H). The minimum in 5T n is achieved when 



5T n = min(2, —J—). 

a n V K 

Let us consider two different regimes. 

1. If < 1 then the minimum is achieved for 5T n > 1 and 

B(k, vq, U,H)>C max { min ( ^(1 - a n ) + ^ + v ) 1 , 

0<h<H/3 [0<a n <l \ H h J J 

where C is a universal constant. Choose h out of the condition Uh/H = n/h, so that 



Then evidently 



B( i; .r l) .U.II) >C {\l^ + vo \ =CvJJ^ + l 



with t c and t u defined in (70). Given ( p.01| ), condition ^ < 1 translates into 



Vq V H T u 

For this computation to apply, we also need to ensure that h < H/6, that is, 

K 

UH £ 1/6 ' 

so that h is in the acceptable range. This is a condition on the Peclet number. All our bounds will 
remain valid if replace the choice of h by Ch with C in some fixed range; we just have to adjust the 
universal constants. Therefore, for presentation purposes we will henceforth require the condition 

m 



Pe = \ > 1. 



2. If rrr > 1 then, since 5T n = min(2, ^\ P^r), we have that 

B(h,vq,U,H)> max min afar,), 
v ' ; ~ o<h<H/6 0<a„<i yv 
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where 

Uh 2 v 1 1 _ \ , 2Ka n , < l / ^ n fe 



Note that 



ma X/ nun f (1 - a n ) + + ^) > C{(Uv qH -^)^ + , ) = C Uo f f ^ ^ + , 

0</i<H/6 0<a n <l y kH h 2 J \\T U J J 

since we can choose h from Uh 2 vo/(K,H) = n/h, which gives 



Uv 



Also 



. (Uh 2 V N , „ [v^k\ ^ „, TT 3 rT _i N i/ 5 ^ /^cA 1/5 

max mm — — ( 1 — a r , 

0<h<H/30<a n <l \ kH 



since we can choose h from Uh 2 vo/(KH) = \JvqkJTi, which gives 

( k 3 H 2 
\U 2 v 



Therefore, 



max min ofa^) > C^nmin < — ) + 1 

0</i<///3 0<a„<i yV n/ - I 1 1 - ' 



Notice that t c /t u > 1 in the regime we consider, since t c /t u is also equal to {k/vqK)^ with h chosen 
according to ( |102| ) or (K/voh) 5 / 2 with h chosen according to ( |103j ). Therefore, we get that the second 
regime is characterized by 

Il = l [^L> l 
t u v ]j H - 

and in this case 

B(k,v q ,U,H) > Cv 



The choice of /i is given by ( |103| ). In this regime, we also need to assume that 

1/5 



so that h is in the acceptable range. This, up to a constant, is a condition 

/ 

H 

which is satisfied provided (|69|) and (p8|) hold. 



< Pe 2 
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Now we are finally ready to give a proof of Theorem 0. 
Proof. Given the parameters k, vq, H, and U, define h according to 



1 kH , , , 

i /k 3 h 2 \ 1/5 

h= i[v^) • T < /T » aL (105> 

These choices ensure that h < H/6 provided that (|68|), (|69|) are satisfied. Consider the smallest 
integer number mi such that (T)j < 04 — (04 — 0i)/lO for all i > 171%. We need to consider several 
options. The first is taken care of by 

Lemma 9 If (T) mi < 04 — (04 — 0i)/5 then the bounds of Theorem hold. Moreover, the same is 
true if for some m, m + 1, we have (T) m , (T) m+ i £ [0i, 64] and 

\(T) m+ i - (T) m \ > (0 4 - #i)/10. 

Proof. Let us consider the case where (T) mi < 04 — (04 — 0i)/5, the other case is similar. If m-i = 2n 
is even, notice that by definition of mi, (T)2 n -i > 04 — (04 — 0i)/lO. Then by the diffusive estimate 
(94) we obtain 



1 „k 



v>^-v n >c~(e 4 ° a2 



3 " - /i 

In the case where t c /t u < 1 putting in h given by (104) gives exactly the bound ( f7ll) of Theorem 0, 
modulo C2^o- But by Proposition ||, we can always add Cvq to the lower bound on V. In the 
case where t c /t u > 1, putting in h given by (|105| ) gives the lower bound C(K 2 [/ 2 i/- 2 t> ) 1/5 = 
Cvo(t c /t u ) 2 / 5 , which is even better than the bound ( |7l| ) (modulo vo). Hence, Theorem ^ also holds. 

If mi = 2n + 1 is odd, then by definition of mi, {T)2n > 04 — (04 — 0i)/lO. Thus temperature 
falls inside the cell C n , and we may use advective estimate of Theorem ||. We can find £oj £1 i n the 
cell C n such that 

3h 3h 



- J k(p)T(p,^)d P = e 4 -^ T ^, - J k( P )T(p,h)dp 



and 



3h 

1 



k( P )T(p,Z)dpe [04 



F J 5 10 

ft 

for every £ between £0 and £1. By the advective estimate (|S^), we get 



V>3V„>C(C./„)^l + ^j (».-».)V 

Direct substitution of the expression ( |104| ) or ( |105| ) for h depending on the value of t c /t u and 
comparison of the above bound with ( |7l"D gives the conclusion of Theorem ^. The proof of the 
second statement of this lemma is parallel to the above argument. □ 

Lemma |9| proves Theorem [?] unless there exists a sequence mi, . . . , (of length at least 8, in 
fact) such that for every mi < m < m-2, 

(T) m € [0i + 94 ~ 91 " 94 



10 ' 10 
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and for m = m\ , . . . , m 2 — 1 



\\T) m +l — {T) m \ < 



10 



Therefore there exists a sequence of at least three (or more) consecutive cells C ni , C n2 such that 
(T) 2ni > 8 A - A(fl 4 _ e x ), (T) 2n2+1 < 0x + A(# 4 _ 9l ), (106) 

and 



(107) 



10 ' 10 J 

for every m such that 2n\ < m < 2n 2 + 1. Assume that there exists a cell C n with m < n < n 2 
where we can find £o so that 



3/i 



k(p)T(p,Zo)dpt [8 



1> C4 



By ( |107| ) and mean value theorem we can also find £i in this cell such that 

3h 

J k( P )T( P , eo ^ e + tf 4 - 

h 

Therefore, the advective estimate (|82|) can be applied in this cell, giving 

Uh 

v oh J v ±y H 

Hence, similarly to the proof of Lemma |9[ Theorem |7] holds in this case. The only case left to 
consider is the case where for every n\ < n < n 2 , for every £ in a cell C n , we have 



V n >C(CJo)(l + ^r) 1 >4-#i) 3 



3/i 



fc(p)r(p,£)dpe [<?i,e 4 ]. 



In this case the second advective estimate (^3|), as well as reactive estimate ( pq ) apply in every cell 
C n such that n\ < n < n 2 . Recall the notation 

(1 — a n )8T n = ST n: a = |(T) 2n+ i — (T) 2n |, a n 6T n = 8T nt d = \{T) 2 (n+i) — {T)2n+i\- 

Following the computation we performed at the beginning of this section we get 

1 Uh 



y>\nyn>cY. 



— y 

vohj H 



l-a n ) + -a n 5T n + v 5T n J ST, 



(108) 



n=n\ n=ni 

Consider the case where t c /t u > 1 (the other case is similar but simpler). Putting the expression 



for h from ( 105 ) into fllOSj ), we get 



ll 2 



v>cJ2 



n=ni 



1/5 



vol—) (1 - (hi) + vo ( — ) al5T n + v 6T n 1 



2/5 



5T n 



Since by f06|) , En= m ST n > 2(# 4 - 0i)/5, we have 



V > Cfn min 

0<a„<l,0<(5T„<2 



1/5 



1 - O n ) + 



2/5 



a 2 JT n + 5T„ 1 



(109) 
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It remains to show that the expression in square brackets in ( |10S| ) is always greater or equal to 
Cvq^Tc/tu) 1 / 5 (we can always add Cvq later to the lower bound for V by Proposition [[]). If a n < 1/2, 
the first term in the sum gives exactly the estimate we need (no matter what is the value of 5T n < 1). 
Hence it remains to consider the case where a n > 1/2. In this case the sum of the second and third 
term in the square brackets is greater than or equal to 

c -((S) 2/5ir » +5r " 1 ) £C "°(S) 1/5 

and this finishes the proof of Theorem [?]. □ 
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